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ABSTRACT 


This investigation was conducted to determine the advantages and disadvantages of using thrust 
vectoring for lateral-directional control and the effects of reducing the tail size of a single-engine 
aircraft. The aerodynamic characteristics of the F-16 aircraft were generated by using the Aerodynamic 
Preliminary Analysis System II panel code. The resulting lateral-directional linear perturbation analysis 
of a modified F-16 aircraft with various tail sizes and yaw vectoring was performed at several speeds 
and altitudes to determine the stability and control trends for the aircraft compared to these trends for a 
baseline aircraft. A study of the paddle-type turning vane thrust vectoring control system as used on the 
National Aeronautics and Space Administration F/A-18 high alpha research vehicle is also presented. 

INTRODUCTION 

Increasing interest in the use of thrust vectoring control systems to enhance or replace aerodynamic 
controls for future fighter aircraft exists. Pitch vectoring nozzles have the potential to increase turn 
performance, reduce drag, and increase short field flight operations through careful nozzle integration 
and optimum vectoring. The benefits of yaw vectoring include directional control power at low speeds, 
at high angles of attack, and during spins and stalls, where the traditional control surface effectiveness 
is substantially reduced. Reduction of vertical tail size by using thrust vectoring control power also 
promises to benefit aircraft observability and to reduce drag and weight substantially; however, the 
addition of the thrust vectoring control system may eliminate some of these benefits. 

Now, several aircraft incorporate thrust vectoring in the control systems. The F-15 short takeoff 
and landing maneuver technology demonstrator and the proposed advanced tactical fighter have pitch 
vectoring nozzles, while the X-31 aircraft and F/A-18 high alpha research vehicle incorporate pitch 
and yaw vectoring through use of postexit vanes. Pratt & Whitney Co. (West Palm Beach, Florida) is 
developing a lightweight, highly efficient, three-dimensional nozzle to replace the current nozzles on 
America’s F-14, F-15, and F-16 fighter aircraft. 


Purpose 

The fighter aircraft of the future will require use of thrust vectoring controls to attain the maxi- 
mum controllability during maneuvers and optimum performance. Thrust-vectoring controls will work 
in conjunction with traditional aerodynamic control surfaces to enhance the longitudinal- and lateral- 
directional flying characteristics of the aircraft. This study describes the significant effects of thrust 
vectoring control systems and reduction of vertical tail size on lateral-directional stability and control 
of a single-engine fighter aircraft. 


Aircraft Description 

For this study, the F-16 single-seat fighter aircraft with the F-100 engine model derivative and four 
postexit thrust vectoring control vanes was employed. The General Dynamics (Fort Worth, Texas) F-16 
aircraft is a small, single-engine, single-seat, multirole, tactical fighter known for its high performance 
and agility (fig. 1). The wing of the F-16 aircraft has a moderate sweep and automatic leading-edge 



flaps and flaperons. The flaperons combine the functions of automatic trailing-edge flaps and ailerons 
and provide pitch and roll control. The all-moving horizontal tail has a small amount of anhedral and 
provides both pitch and roll control through symmetric and differential deflections. The vertical tail 
and two ventral fins provide directional stability, while rudder deflections on the vertical tail provide 
directional control (refs. 1 and 2). 

The F-16 aircraft incorporates an analog flight control system with a force- sensing sidestick con- 
troller and associated feedback loops, gain schedules and limiters for angle of attack, g loading, roll and 
yaw rates, and rudder pedal force. The current roll and yaw rate limiters automatically engage at high 
angles of attack and are coupled to the roll and yaw rates and the angle of attack (ref. 3). 

The engine of the F-16 aircraft used in this study is a modified version of the currently deployed 
F-100 engine known as the F-100 engine model derivative and was used in the F-15 highly integrated 
digital electronic control aircraft at the NASA Dryden Flight Research Facility. The F-100 engine 
model derivative (company designation PW 1 128/Pratt & Whitney, West Palm Beach, Florida) is a low 
bypass, twin-spool afterburning, turbofan engine (ref. 4). The primary nozzle throat area is controlled 
by digital electronic engine control. Aircraft trends indicate that the thrust vectoring and reversing 
aircraft will employ digital electronic control systems such as highly integrated digital electronic control 
and full authority digital engine control and will allow an improved engine performance, a larger flight 
envelope, a faster stall-free engine acceleration, a lower specific fuel consumption, and a safer operation 
as it reduces the pilot’s workload during takeoffs and landings (ref. 5). 

For this report, an engine simulation program that predicts the engine performance through com- 
ponent characteristics was used. The F-100 specification model is a steady-state, aerodynamic, math- 
ematical model of the F-100 engine model derivative turbofan engine used to obtain the predicted net 
thrust, gross thrust, and nozzle pressure ratios at specified altitudes and Mach numbers (ref. 4). 

A thorough analysis of the thrust vectoring control system being installed on the F/A-18 high alpha 
research vehicle at NASA Dryden was used to build the thrust vectoring control system model used in 
this study. The three nickel alloy Inconel (International Nickel Company, Huntington, West Virginia) 
postexit vanes that were cold-jet tested at the NASA Langley Research Center (Langley) for the F/A-18 
high alpha research vehicle were used to model a four-vane system proposed for the F-16 aircraft. 

Project Scope 

The lateral-directional linear perturbation analysis of the F-16 aircraft with various tail sizes was 
performed at a wide range of altitudes and speeds to determine the stability and control envelope. 
The aerodynamic data were obtained by constructing several models of the F-16 aircraft using the 
aerodynamic preliminary analysis system panel code. Postexit vane efficiency and effectiveness were 
obtained by analyzing the thrust vectoring control system cold-jet tests of the F/A-18 high alpha research 
vehicle at NASA Langley. Next, a linear perturbation analysis program called Stability 1 was constructed 
to obtain the stability and control characteristics. This program incorporated the cold-jet efficiency and 
effectiveness coirelations, the Aerodynamic Preliminary Analysis System II aerodynamic coefficients 
completed for this study, and the results of the F-100 engine specification deck and NASA Langley’s 
wind-tunnel tests of thrust induced effects on the high alpha research vehicle. 
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Organization 


This report is comprised of sections dealing with the advantages of thrust vectoring control systems, 
a discussion of thrust vectoring control systems, the aerodynamics of the aircraft, the stability analysis 
results, and the conclusions and recommendations. 

NOMENCLATURE 


Abbreviations 

Aerodynamic Preliminary Analysis System II 

advanced tactical fighter 

digital electronic engine control 

engine model derivative 

full authority electronic control 

high alpha research vehicle 

highly integrated electronic control 

lateral coefficient spin parameter 

leading-edge flap 

National Aeronautics and Space Administration 

nozzle pressure ratio 

postexit vanes 

poststall maneuver 

F-100 engine specification deck 

Stability- 1 computer program 

short takeoff and landing 

takeoff gross weight 

thrust vectoring control system 
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ADVANTAGES OF THRUST VECTORING CONTROL SYSTEMS 

The future air warfare arena will be quite demanding and will require improvements in agility, 
control, performance, short takeoff and landing (STOL) capabilities, and stealth. The use of thrust 
vectoring is being considered as a way to accommodate these demands. 


Agility 

Aircraft agility can be defined as the aircraft’s maneuver and control capabilities, that is, the ability 
to maneuver and point quickly, to obtain first-shot advantage over an opponent, to sustain high-tum rates, 
accelerations, and decelerations, and to perform offensive and defensive maneuvers at will (refs. 6 and 
7). Agility is currently measured by pitch rate, load factor, turn performance, roll rate, level yaw rates, 
and acceleration and deceleration (ref. 8). 

History shows that aircraft designed during times of war tend to have greater agility than those 
designed during times of peace (ref. 9). The F-16 aircraft was designed for high maneuverability; 
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however, its wing loading has increased substantially as avionics and weapons systems have been added 
over the years. To compensate for the loss of maneuverability, a modified F-16 aircraft, the Agile Falcon, 
with an increase in wing area and the addition of thrust vectoring has been proposed. Maneuverability is 
important; however, improvements in maneuverability must be accompanied by increased controllability 
or overall agility will suffer. During the Korean conflict, the MiG 15 aircraft had greater maneuverability 
than its adversary, the F-86 aircraft; however, the F-86 aircraft had greater controllability and was able 
to rack up an impressive ratio of air victories (ref. 8). 

Controllability is the quickness of ease to change the flightpath of the aircraft from moments and 
forces caused by the control surfaces. In addition, controllability gives the option of stepping out of the 
limits of optimum efficiency or to an arena of stability not considered normal. Unfortunately, current 
flight controls of the F-16 aircraft try to prohibit unstable high angle-of-attack maneuvers through use 
of pitch, roll, and yaw rate limiters. The goal of increasing controllability is to increase the envelope 
of positive control for the pilot, especially at high angles of attack or sideslip. A solution for increased 
controllability demands is for a thrust vectoring control system (TVCS) to provide the required forces 
and moments (ref. 6). 

The F-15 short takeoff and landing (STOL) maneuver technology demonstrator (MTD), which 
incorporates a two-dimensional pitch TVCS, is designed to have increased turn performance, 33 percent 
improvement in pitch rate and 53 percent improvement in roll rate because of the TVCS (ref. 10). As 
pitch and roll are enhanced, yaw rate promises to be substantially increased because of yaw vectoring, 
especially at high angles of attack where roll and yaw rates are seriously degraded due to the control 
surfaces being blanked out by the disturbed flow of the aircraft wing. 

Poststall Maneuvering 

There has been a recent trend toward the exploration of poststall maneuvers (PSM) through studies 
and flight aircraft. The Russian demonstration at the 1989 Paris airshow of a production Su-27 perform- 
ing a Pugachev’s Cobra, where the aircraft obtains an angle of attack greater than 90° without loss of 
control, has reemphasized the importance of exploring this regime (ref. 11). Currently, NASA Dryden 
Flight Research Facility (Dryden) is flight testing the F/A-18 high alpha research vehicle (HARV), which 
has achieved studies in angles of attack of up to 55°, and is dedicating flights of the X-29 aircraft to 
the study of forward-swept wings at high angle of attack (ref. 12). 

Poststall maneuvering increases close-combat effectiveness by allowing the advantage of quicker 
entry into and longer maintenance of firing position, more first-shot opportunities, and fewer counter 
hits. Poststall maneuvering studies indicate that it is easier and quicker to switch targets, deploy missiles, 
and dictate tactics throughout the entire speed regime. These advantages lead to an exchange ratio of 
2:1 in dual-combat studies against opponents of equal conventional aircraft (ref. 7). 

The Herbst maneuver is traditionally associated with a short excursion into the poststall regime and 
then a return to normal flight. As a combat technique, the Herbst maneuver is characterized by its high 
angle of attack; its short duration, an average of 5 sec; and its use in 10 percent of engagement time, 
lower overall g levels, and lower overall speeds during air-combat simulations (ref. 7). Longer poststall 
maneuvers include a helicopter maneuver that allows pilots to track the target by yawing instead of using 
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a coordinated turn and a fully controlled tail slide maneuver that allows pilots to track an opponent that is 
attacking from above. Poststall maneuvers traditionally have high pitch rates and include simultaneous 
rotations in roll and yaw (ref. 13). 

H. Ashley found that in real-time simulation studies that high-angle-of-attack requirements are 
highly oscillatory and unsteady and require large forces and moments over a wide range of attitudes 
(ref. 14). He suggests that ±10° of total jet deflection of thrust are needed to supply control because of 
the inadequacy of the aerocontrols. 

For a pilot in an offensive position, pitch and yaw vectoring can be a powerful maneuver depending 
on the type of opponent, the opponent’s flight regime, and the ability to optimize maneuvers. If both 
aircraft are not in a low-speed flight regime, entering a poststall maneuver, where thrust vectoring would 
be employed while the opponent maintained a higher energy state, would be a great disadvantage. 
Furthermore, if both aircraft are in low-speed flight, the pilot must perform an optimum maneuver to 
gain the maximum advantage. Pitch and yaw vectoring tends to be less effective at the lower altitudes 
because the advantages of turn rate are negligible. The advantages in turn rate and agility zone increase 
as altitude increases. It is advisable to incorporate yaw vectoring into the controls and to have the 
option of turning on or off the pitch vectoring at the pilot’s discretion. Simulation studies showed that 
pilots require larger positive acceleration because of gravity (g ) loading and upward pitch vectoring 
than downward pitch vectoring. Because poststall maneuvers and short takeoffs require large vector-up 
deflections, while the downward vector requirements for optimum performance are usually less than 
10°, tilt the nozzle to increase the upper limits of the vectoring envelope. 

Stealth 

The one area in which great benefits can be achieved through the use of thrust vectoring is that 
of stealth. Because radar cross-section is a measure of the observability of the aircraft and is an 
logarithmic function of the profile area presented by the aircraft, small reductions in this area realize 
great benefits. From a side view of the aircraft, the fuselage, cockpit, and vertical stabilizer are the 
main contributors of observability (ref. 15). The current stealth modifications to the F-16 aircraft include 
radar absorbing materials in the intake and a reflective gold-covered canopy. Because the fuselage size 
is relatively constant, it can be concluded that the reduction of the vertical stabilizer has the largest 
room for improvement (ref. 16). The addition of a TVCS may nominally increase the observability 
of the aircraft, but a large reflective vertical tail is harder to mask than a faired TVCS. Quantitative 
improvements in radar observability are very difficult to predict and require thorough testing. The 
amount that the tail size can be reduced is largely a consideration of the dynamic stability and control 
of the aircraft. 


Thrust Vectoring Concepts 

In the past, most TVCSs have been developed for pitch vectoring concepts such as the F-15 S/MTD 
in a trend toward developing short field capability technology and performance optimization. In recent 
years, the interest in PSM has increased the investigations of pitch and yaw vectoring concepts. Similar 
in application, pure yaw vectoring concepts are less beneficial in overall aircraft performance than pure 
pitch vectoring two-dimensional nozzles. 
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There are two types of nozzle classifications: the nonaxisymetric nozzle and the conventional ax- 
isymmetric nozzle. The nonaxisymmetric nozzle concepts typically employ sidewall flaps or rectangular 
ports to affect the direction, while the typical axisymmetric nozzles mechanically turn all or part of the 
nozzle or employ postexit vanes. With nonaxisymmetric nozzles, rounded or rectangular nozzles are 
simpler to implement vectoring and have fewer moving parts than an axisymmetric nozzle with variable 
throat control. At low Mach numbers, twin-engine aircraft with nonaxisymmetric nozzles have lower 
drags, but as the Mach number increases to the transonic region, the axisymmetric nozzles exhibit 
lower overall drags. In addition, nonaxisymmetric nozzles may be more sensitive to adverse tail effects 
than axisymmetic nozzles at transonic and supersonic Mach numbers where these effects are significant 
(ref. 17). 

The TVCS have several disadvantages, the primary of which is the reduction of thrust performance 
and the lost nozzle efficiencies. Losses occur in nozzles because of skin friction, internal flow separation, 
exhaust flow divergence, under- and overexpansion, and temperature and thrust distortions. Turning the 
gross thrust from axial direction also causes losses in nozzles. For postexit vanes, additional skin friction 
and pressure losses are caused by the thrust vectoring hardware (ref. 18). 

Two-dimensional, rectangular, nonaxisymmetric nozzles integrate well with tapering nearly rectan- 
gular aftbody fuselage contours, thereby offering a potential reduction in the nozzle airframe drag during 
cruise (ref. 18). For twin-engine aircraft, this configuration would have favorable effects with nozzle 
airframe-induced drag. With careful attention, nozzle boat tailing can be eliminated in the immediate 
vicinity of the nozzle interfairings, tail-support boom structure, and empennage surfaces, and, therefore, 
in minimizing nozzle aft end-flow separations. Through the positive effects of optimum vectoring, drag 
reductions of up to 13 percent have been obtained (ref. 19). Drag reductions from axisymmetric and 
nonaxisymmetric nozzles can be made by careful attention to the wing-body geometry and its interac- 
tion with the flowfield. Success depends on the local contouring and the rate of area change of aftbody 
nozzle area. The drag reductions are the most significant by use of the thrust-induced, supercirculation 
effects for enhancing effects at maneuver conditions. The very dramatic changes between the jet-on and 
jet-off conditions, even when unvectored, display complex characteristics that will require greater levels 
of understanding than conventionally employed in nozzle design and airframe integration (ref. 20). 

The structural efficiency of nonaxisymmetric nozzles is less than that of axisymmetric nozzles. 
Further weight gains can be attributed to the fact that the engine case structure must be strengthened 
to withstand the vectoring. Takeoff gross weight penalties from 5 to 7 percent occur because of the 
structural effects and cooling system weights (ref. 18). In one study, the weight of a prototype three- 
dimensional thrust vectoring nozzle on an RB199 engine was 254 lb (ref. 13). For the postexit vanes of 
the F/A-18 HARV, the total weight gain was 1200 lb for the TVCS for a 3.33 percent increase in takeoff 
gross weight (TOGW). In addition, there was an 800 lb increase or 2.22 percent greater TOGW for 
the counterbalance weight to maintain aircraft center of gravity (c.g.) position. For the postexit vanes 
(PEV), the TVCS comes to 200 lb/vane. 

At this time, Pratt & Whitney is developing an axisymmetric balanced beam pitch and yaw vectoring 
nozzle for installation on fighter aircraft with very little weight increase. This new design incorporates 
many standard, nonvectoring, engine nozzle components and can have ratios of thrust deflection over 
nozzle deflection greater than unity. When this system is operational, the PEV thrust-vectoring system 
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will be obsolete. For this study, postexit vanes are used because of the availability of data and to 
demonstrate the maximum effects of implementing a thrust-vectoring system. 

THRUST VECTORING CONTROL SYSTEM 

The choice of an axisymmetric nozzle with paddle-type postexit vanes was based on modifying 
the F-16 aircraft using the current technology of the F/A-18 HARV and X-31 aircraft programs. These 
projects chose to implement PEV because installation is relatively simple, and it is time and cost 
effective in demonstrating the effects of thrust-vectoring systems. The PEV are separate from the 
nozzle, so a conventional axisymmetric nozzle could be used without reducing nozzle performance or 
adding complicated, costly vane cooling. 


Postexit Vanes 

From a study by B. Berrier and M. Mason, it was found that the postexit vanes were competitive 
with other multiaxis, nonaxisymmetric systems even though large losses were incurred (ref. 21). In 
addition, spoon-shaped vanes were more effective than single-curvature vanes. On the other hand, 
single-curvature vanes were more efficient than flat vanes. The type of nozzle used is of little effect 
on the postexit vane performance; therefore, for this study, the standard F-16 aircraft nozzle can be 
used. The research also showed that using four vanes rather than three provides a larger thrust envelope 
(fig. 2). Berrier and Mason’s study also confirmed the results for the F/A-18 HARV that increasing vane 
area relative to the nozzle exit area increases the resultant thrust vector angle (fig. 3). 

The F/A-18 project has conducted static performance testing of the postexit vanes at the NASA 
Langley 14 x 22 ft subsonic tunnel (ref. 22). The three- vane configuration was chosen to conform with 
the aftbody structure and resulted in two differently sized vanes at three different relative angles (figs. 4 
and 5). 

The vanes were tested at both maximum afterburner and military power and were positioned at the 
following angles: —10, —5, 0, 5, 10, 15, 20, 25, and 30 deg; however, the F/A-18 HARV vane-deflection 
angle is structurally limited to 25°. The test results were compiled and used to draw correlations for 
the size and position of the thrust vectoring system to be used in this study (appendix B). 

The TVCS vane for the F-16 aircraft in this study was chosen to use highly effective, double- 
curvature vanes, the same as used on the F/A-18 HARV. The F-16 aircraft will employ a convergent- 
divergent standard nozzle. In addition, the four vanes will be placed symmetrically on the aftbody for 
a larger control envelope and simpler control system than a three-vane configuration permits (fig. 2). 
From the results of the F/A-18 HARV static tests, correlations between the jet deflections, efficiency, and 
direction with thrust vectoring control deflection were made to predict the performance characteristics 
of the F-16 aircraft. 

Figures 3 and 6 show that for single-vane deflections the larger the vane area, the greater the 
thrust deflection. In addition, the maximum afterburner power setting is more effective than the military 
power setting for vane deflections between 10 and 25° (figs. 3 and 6). Nozzle pressure ratios significantly 
increase the thrust deflections, especially at smaller vane-deflection angles. This increase occurs because 
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of the expanded flow contacting the vanes at smaller deflections (fig. 7). The maximum NPR effects 
occur from 15 to 25° for military power and 10 to 15° for maximum afterburner with NPRs of 5 and 
6 having the greatest differences. 

The two-vane deflections have more vane area relative to the nozzle exit area and, therefore, have 
greater thrust deflections compared to single vanes. It was found, however, that significant interference 
effects occurred in large deflections and were highly dependent on power settings, vane sizes, positions, 
and deflections. Figure 8 shows that the A and C vane configurations, which are positioned 138.5° apart, 
resulted in less thrust deflection than the smaller B and C vane area combination, showing the importance 
of relative angle to vane interference. In figures 9 and 10, it is seen that vane interference increased with 
nozzle pressure ratio and that all vane configurations were affected. Comparison of figures 8 and 10 
with figures 11 and 12 reveals that the interference effects for maximum afterburner power are greater 
than the military power setting and effectiveness decreases further with nozzle pressure ratio. The two 
vanes deflected the jet far enough to impinge the jet on the inactive vane therefore eliminating further 
jet deflection and incurring high losses. At higher nozzle pressure ratios with the flow expanding as it 
leaves the nozzle, the jet-deflection performance was further reduced. Figure 12 shows that vane size 
and position have a major effect on jet deflection and interference. Note that inactive vane position, 
size, and deflection played a role in interference. 

The proposed F-16 aircraft would use the standard HARV TVCS because of the availability of 
the data. Future designs should incorporate asymmetric vane deflections for opposite vanes or position 
the inactive vane farther from the nozzle so as not to impinge on the vectored engine thrust Either 
configuration would reduce interference and increase jet-off effectiveness. 

The direction of the vector for the vanes reacted as expected. When a vane at a specified radial 
position is deflected, the thrust vector radial direction should be opposite the vane radial position. For 
very small values of total vane deflection, the readings produced very small vector deflections and 
showed considerable error. On the other hand, as the vane deflection increased to 15°, where deflection 
becomes significant, the radial directions approached the expected angle (fig. 13). Similar correlations 
can be found in the case of two vanes deflecting at once. The direction is a function of vane radial 
position and radial vane size. When both vanes have the same deflection angle, for example, 

ARC\ flj ■+■ ARC 2 Q 2 
= ARC, + ARC 2 

where 


ARC = radial vane size 

fi = radial vane position angle 


This equation assumes that similar vanes are used, and there is no third-vane interference. 

The efficiency of the postexit vanes is the total magnitude of the deflected thrust vector over 
the magnitude of the undeflected thrust vector. Efficiency for TVCS vanes is poor when the vanes 
are deflected beyond 20° (figs. 14 and 15). In figures 16 and 17, the smaller area B and C vane 
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combination is slightly more effective than the other combinations, and there is a crossover of efficiency 
with nozzle pressure ratios because of inactive vane interference. The nozzle pressure ratio decreases 
efficiency slightly, and the maximum afterburner case is substantially more efficient than the military 
power setting. 

The efficiency and deflection of the postexit vane thrust vectoring system is highly dependent 
on active and inactive areas, placement, vane deflections, and nozzle pressure ratios. The maximum 
afterburner configuration has the greater deflections, efficiencies, and interference effects compared to 
the military power setting. 


Thrust Vectoring Simulation 

The thrust vectoring simulation subroutines used to determine the deflection, direction, and effi- 
ciency are based on the thrust vectoring data from the F/A-18 HARV cold-jet tests. The single-vane 
thrust deflections curves with respect to the vane deflections were approximated using fourth-order poly- 
nomials at each nozzle pressure ratio for both military power and maximum afterburner. The efficiency 
correlations with respect to the single vane deflections used third-order approximations. The two vane, 
equal deflection results are modified single vane deflection models which account for the interference 
of two differently sized vanes at a relative radial angle. The vane-deflection and efficiency-correlation 
programs are reasonable approximations within the angles, deflections, and sizes of the F/A-18 HARV 
configuration. Better results for single-vane deflections are achieved because the correlation programs 
were directly interpolated from the single-vane data, and no interference effects were present. 

The F-16 radial vane width of 60° was based on the F/A-18 HARV ratio of vane radial width to 
that of the total radial degrees. Vane relative angles of 90° were chosen to simplify control laws, obtain 
maximum effectiveness and efficiency, allow vane placement without restricting horizontal tail and 
rudder deflections, and gain maximum thrust vectoring in pitch or yaw with two equal vane deflections. 

Aircraft Configuration 

The F-16 fighter used for the study incorporated an F- 100 EMD engine with DEEC control. Four 
Inconel, double curvature, postexit vanes were positioned at ±45° and at ±135° for pure pitch or 
yaw vectoring with two vane deflections (fig. 4). The vane control power was derived from linear 
interpolation of two equal-vane deflections of 20° assuming that the vane- zero position was at 0° vane 
deflection. The vane deflection, direction, and efficiency models are based on the F/A-18 HARV and 
assume no differential vane deflections and no control surface limiters. The leading-edge-flap schedule 
for the F-16 aircraft was used, but the horizontal tail was untrimmed and at 0° deflection. The c.g. 
positions are maintained at 30 percent of the mean aerodynamic chord for all configurations. The 
weights and inertias are as follows: 
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Wi 0 = 16,3301b 
Ixxq = 9,470 slug/ft 2 
I yyo = 49,450 slug/ft 2 



Izz Q = 56,760 slug/ft 2 
Ixzq = 310 slug/ft 2 

From Nicolai (ref. 23) the weight of a vertical stabilizer for a conventional aircraft with subsonic and 
supersonic capabilities is 

W tvl = 0.19 7 1014 

7= (l + ^° 5 (W to n) 0363 (S n ) hm M maXa l _ om L vt 
(l + J^) 0 ’ 217 AR n a337 (l + A vt ) a363 (cos A^)- 0484 


With the following results, the weights for the vertical tail and two ventral fins are 

W tvl = 3001b 


^tventral fins 


= 1061b 


where = 0 for a fuselage-mounted horizontal tail. The weights for the proposed F-16 TVCS and 
counterbalance are W^cs = 600 lb and Wt cb = 600 lb. The inertias are defined as follows: 

I XX = J (y^ *F 2 2 )6?7l 

Iyy = J (x 2 + Z 2 )6m 
hz = J( x2 + y 2 )6m 
Ixz = [ (xz)6m 


J v ' 

Ixy ~ lyz — 0 


hytv + cb 

* zz tv C J, 

Ixx yt 100% 
Iyy vi ioo% 
I zz vl 100% 
Ixz vl 100% 


= 14,906 slug/ft 2 
= 14,906 slug/ft 2 
= 84 slug/ft 2 
= 3, 143 slug/ft 2 
= 2, 387 slug/ft 2 
= 559 slug/ft 2 
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The reduction of tail height was assumed to be proportional to the reduction in the vertical tail weight. 


AERODYNAMIC PRELIMINARY ANALYSIS SYSTEM 

To obtain the coefficients for the F-16 aircraft flight envelope, the Aerodynamic Preliminary Anal- 
ysis System II (APAS) was used. The APAS is a computer code on the NASA timesharing network that 
uses a unified distributed panel method and incorporates potential flow theory with edge considerations 
at subsonic and supersonic speeds in the analysis of three-dimensional configurations having multiple 
nonplanar surfaces of arbitrary platform and of bodies of noncircular contour. These surfaces are ana- 
lyzed statically, rotary, and longitudinally. Then, longitudinal and lateral-directional characteristics are 
derived. In preliminary design, the finite-element analysis combined with realistic assessments of its 
limitations, is a valuable tool (ref. 24). 

The jet-flap and jet-interference effects options on APAS were not used because of inoperability. 
As a result, the F-16 aircraft intake was faired in for the APAS configuration (fig. 18), and the effects of 
unvectored engine thrust were assumed to be small. Singularities in the system occurred when operating 
at Mach numbers 1.0 and 1.2. As a result, the aerodynamic data were taken at Mach numbers 0.2, 0.4, 
0.6, 0.8, 0.9, 1.1, 1.25, 1.4, 1.6, and 1.8. For this study, the c.g. position was assumed to be at 30 percent 
of the mean aerodynamic chord at the frontmost position of the c.g. envelope for the F-16 aircraft. 

Stability Coefficients 

One objective of this study was to examine the effects of reducing the vertical tail size on aircraft 
stability and control. The five configurations that were studied with the APAS code are presented 
in figure 18. The APAS studies produce stability coefficients that followed the expected trends with 
respect to Mach number for all but low supersonic speeds (refs. 25 and 26). The variations of the yawing 
moments, rolling moments, side forces with sideslip angle, and roll and yaw rates showed a decrease 
in the magnitude of the coefficients as tail size was reduced (figs. 19 to 30). The large increases in the 
coefficients at Mach numbers 1.1 through 1.4 were much larger than expected and are considered to be 
in error. 

As was expected, directional stability decreased as the tail size was reduced. Directional stability 
and rudder control power were reduced by 35 percent with configuration B and 70 percent for configu- 
ration C (fig. 19). The aircraft appears to have become statically unstable, Cnp < 0.0, between 20 and 
30 percent tail height. Although the magnitudes of the lateral-directional stability coefficients do not 
match well with those presented in a study of the prototype F-16 and the F-16 aircraft with the c.g. at 
35 percent (refs. 26 and 27), the ratios of the coefficients must be accurate for the linear perturbation 
analysis to be correct 

The longitudinal coefficients for the baseline aircraft were very similar to those expected; but, when 
the tail size was reduced, some uncharacteristic variations occurred. The longitudinal coefficients should 
have remained the same; however, in the lateral-directional analysis of APAS, a small longitudinal 
variation error occurs (figs. 31 and 32). Because the longitudinal characteristics are of secondary 
concern for the small perturbation analysis and are only used to obtain the lift coefficient, the results 
are considered acceptable. For the stability analysis of the F-16 aircraft, the 100 percent vertical tail 
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configuration was used for all longitudinal coefficients. Comparison of lift slope coefficients between 
APAS at 30 percent c.g. and wind and water tunnel models at 35 percent c.g. reveals an average 
difference of 8 percent (ref. 28). 


Vertical Tail Considerations 

The vertical tail provides directional stability and control and must be designed to accommodate 
safe operations such as high- and low-speed flight, takeoffs and landings in crosswinds, one-engine-out 
stability, high angle-of-attack flight, and recovery from spins. Supersonic directional stability is a result 
of small differences in the tail surfaces and fuselage stability characteristics. For high-speed flight, the 
fuselage stability characteristics do not change; however, the vertical tail surfaces may have decreased 
effectiveness because of elastic deformations and become critical at higher Mach numbers (ref. 29). 
This problem is accentuated at high angles of attack and at both subsonic and supersonic flight. The 
F-16 aircraft has an all-moving horizontal tail and two ventral fins; however, the APAS and stability 
analysis codes used in this study do not account for the aeroelastic effects caused by high-speed flight. 

Since supersonic planes fly at high altitudes, the stability derivatives become increasingly ineffective 

as air density decreases, and the inertial effects become more important. The stability of the aircraft 

C n C n 

depends on the parameters -J-&- and T^~- Since Cnp decreases with increasing Mach number, larger 
tail areas, or an artificial substitute are required. Because the postexit vane thrust vectoring system 
adds more mass and inertia to the aircraft than the reduction in tail size, these effects become a further 
decrease in high-speed directional stability. 

TVCS Aerodynamic-Induced Results 

For axisymmetric and nonaxisymmetric all-moving nozzles, there was found to be no aerodynamic 
cross coupling of control forces or moments (ref. 30). For postexit vanes, however, at low vane and low 
jet-deflection angles, there were large discrepancies in the jet direction because of the very small values 
of vane forces (fig. 12). Interference effects between the vanes and vane structural bending are other 
factors which may have caused small amounts of error in jet direction. For the model employed for 
this study, it was assumed that there was no cross coupling between the thrust vectoring control vanes. 

From the studies done for the F/A-18 HARV project conducted at the NASA Langley 30 x 60 ft 
wind tunnel on a 16 percent model of the F/A-18, the induced aerodynamic effects of the TVCS were 
obtained (ref. 22). Drag and lift were decreased by vectoring up and increased by vectoring down 
as expected. These changes are analogous to a blown flap, increasing lift with blowing down and, 
thereby, increasing induced drag. Similar results for lateral-directional vectoring were observed. The 
TVCS-induced effects for the F-16 aircraft have values that are less than 4 percent of the magnitude of 
the yawing moment and side force stability coefficients at Mach 0.2 and at 15° angle of attack. The 
correlations between wing area, wing span, number of engines, and number and size of vanes were used 
between the F-16 and F/A-18 aircraft (appendix B). 

Jet-off effects of postexit vane deflections are significant for differential vane deflections. Although 
these effects are not incorporated in the F/A-18 HARV or in this study, such jet-off effects should be 
considered for future thrust vectoring aircraft research. The postexit vanes act as aerodynamic surfaces 
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such as rudders, elevators, or speed brakes. The differential deflections of thrust vectoring vanes will 
eliminate losses caused by interference and increase the control power available since the vane area can 
be 25 percent of that for the rudder. 


Drag Effects 

The largest drag penalties at transonic speeds result from the effects of the horizontal and vertical 
tail surfaces. On studies of the F-15 S/MTD model, removing the vertical and horizontal tails reduced 
the overall drag from 15 to 16 percent, while vertical tail removal alone resulted in a 10 percent drag 
reduction in the transonic regime. Slightly more than half of the drag reduction was because of the skin 
friction drag, and the rest was because of the adverse interaction with the aftbody (ref. 17). 

As expected for the APAS results, the aircraft drag reduced proportionally as tail size was reduced. 
The largest drag reductions for APAS were made in the skin friction drag (fig. 33). A 14.5 percent 
reduction in skin friction drag occurs in the removal of the vertical tail and ventrals. An 11.3 percent 
reduction occurs with removal of the vertical tail. The wave, base, and vortex drags are relatively 
unaffected by the change in tail size. Total drag values at Mach 1.1 and Mach 1.25 appear to be much 
larger than wind tunnel results and are considered to be in error. 

The drag coefficients were calculated at 0° angle of attack with no surface deflections and, therefore, 
do not include the trim drag, the drag caused by control surfaces, or the unvectored jet-induced drag. 
Skin friction drag can be reduced by removing the vertical tail. Since skin friction drag accounts 
for a little more than half the drag caused by the tail, higher drag savings could be expected. By 
removing the vertical tail, substantial drag reduction could be achieved if drag as a result of the TVCS 
is not significant. Drag increase because of the postexit vanes hardware for 0° angle of attack and 0° 
horizontal tail deflection is calculated to be 0.0027. For the F-16 aircraft, this increase is 13 percent of 
the minimum drag at Mach 0.2 and 6 percent at Mach 1.8. 

The drag penalties of the postexit vanes control system, benefits of the reduction of tail size, and 
optimal pitch vectoring for a trimmed aircraft must be compared to determine the performance gains 
of the aircraft In this case, because the APAS drag results are not considered accurate enough for a 
proper comparison, these calculations were not performed. 

Spinning 

There are two areas in which thrust vectoring could help eliminate the danger of spins: spin 
resistance and control power. Spin resistance is the susceptibility to spiral with and without rolling 
inputs and has significant impact on the handling of the aircraft (ref. 23). The following equation 
describes the region of spin susceptibility for nonrolling turning maneuvers without lateral-directional 
inputs: 

Adynamic = Cn P ~ 
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A negative C na . indicates spin susceptibility. For accessing the spin susceptibility in turning 

^dynamic ** 

maneuvers where lateral-control inputs are introduced, the following lateral coefficient spin parameter 
(LCSP) is employed: 

LCSP = 

For a more detailed discussion of LCSP, see reference 23. When the vertical tail size is reduced, at 
high angles of attack or at high Mach numbers, directional stability of the aircraft, C n/3 , decreases; 
therefore, susceptibility in both the parameters increases. Thrust vectoring control for yaw will not be 
beneficial in improving the LCSP; however, a twin-engine aircraft with thrust vectoring for roll control 
will improve the resistance parameter for rolling maneuvers. A suggested modification for the LCSP 
parameter could be 

LCSP = ^-0^+%^) 

Since thrust vectoring roll would not induce any yawing moment unless there was a relationship similar 
to the rudder and aileron interconnect, yawing moment caused by thrust vectoring would have no effect 
on spin stability. For a single-engine aircraft, such as the F-16, this multiengine characteristic would not 
be a benefit. Relaxing the stability of the aircraft by decreasing the vertical tail size to approximately 
25 percent of its original height will decrease the spin resistance parameters to zero. These resistance 
parameters represent the spin susceptibility of an aircraft without the additional stability resulting from 
control systems with rate and acceleration feedbacks. 

The rudder is the traditional control surface to be employed in a spin to stop aircraft rotation; 
however, careful attention must be paid to the design and placement of the horizontal and vertical tail 
surfaces to ensure that the rudder will be effective. During a spin, the vertical tail surface can be made 
ineffective by the disturbed flow from the wing and horizontal tail. Since thrust vectoring depends on 
engine settings and not on the aftbody flow field of the aircraft, the potential thrust vectoring control 
power will be nearly constant. In addition, the pilot may have the ability to fly in this regime with such 
control as to never enter an unintended spinning state. 

Angle- of- Attack Considerations 

For the F-16 aircraft, static lateral-directional stability remains positive for all angles of attack 
within the limiter boundary of 29°. Directional stability was high below an angle of attack of 20°, 
decreased to one-fourth of its initial value at 26°, and became neutral at an angle of attack of 28° 
(ref. 2). For this study, the angle of attack was never assumed greater than 20° for steady-level flight 
(fig. 34). The importance of yaw and pitch vectoring power with increasing angle of attack is that the 
required control power grows exponentially for a given roll rate (ref. 8). 

The change in angle of attack requires that the values for inertia be transformed to the correct axis 
(ref. 31). For example, 
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The changes that occur in the inertias with the angle-of-attack variation are seen in figure 34. 

Comparisons of Vertical Tail Sizes to TVCS 

The most accurate way to compare the control effectiveness of the thrust vectoring with the vertical 
tail is by calculating the yaw force and yawing moment coefficients. In this case, total net thrust of 
the airplane is defined as T ne t = Tgr 0SS — D^m cos a. For Mach numbers less than 0.4 to 0.5, the ram 
drag is small leaving T ne t ~ Tgross- Since the amount of thrust vectoring available to the pilot during 
steady-state flight is limited to the net thrust output by the engine and equal to the drag of the aircraft, 
an appropriate measure of thrust vectoring power can be derived assuming that the Mach number < 0.4 
to 0.5 and at small angles of attack (ref. 15). If D — Tcosa, then T = Deriving the yaw force 
and yawing moment yields 


*TVCS = Tnet sin 8 tv = 


likewise, 

N TVCS = Tnet ^TVCS sin S tv = a" ^ 


Converting the yaw force and yawing moments into coefficients gives 

n _ YTvrs _ T nM sin fry _ D sin 
°I/tvcs - qsb - qSb “ qSbcosQ 

Similarly, 


n — Ntvcs _ Tnetl-rvcs sin 6 n _ D^tvcs sin i$tv 
° n TVCS - qsb “ qsb “ qsbcosa 

Converting the drag and thrust to coefficient form yields 

n — lives _ Cn sin <$t v 
°j/TVCS qSb=Gr sin 6tv cos a 

likewise, 

Cn„cs = = Curves Siu6 tv = SDh 


At high angles of attack and at high Mach numbers, the thrust vectoring, yawing moment, and force 
coefficient equations must be defined as functions of the gross thrust. The yawing force and moment 
coefficients in terms of gross thrust are derived as follows: 


Vxvcs = Tgross sin 6t v 
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likewise. 


NtVCS = 2gross^TVCS s i° f>tv 


and. 


<Ws = ^ = = C Ttm , sin 

Similarly, 

C„„cs = ^ = = C Tvm J tvcs sin 

For this study, the flight control power was determined for three cases: maximum afterburner, 
military power, and steady level flight The steady level flight control power was determined for all 
Mach numbers and altitudes by calculating aircraft drag and then interpolating for the gross thrust from 
the F-100 EMD data assuming the net thrust equalled the drag. 

The control power of postexit vanes compared to that of the rudder control power is presented in 
figures 35 and 36. In figure 35, the sea level yawing moment available with thrust vectoring is very 
small compared with rudder for the majority of the speed regime except at Mach 0.2, where thrust 
vectoring power is nearly double that of rudder. The yawing moments at 40,000 ft are much smaller 
than at sea level, but thrust vectoring constitutes a much larger proportion of the total control power 
(fig. 36). The magnitudes of the rudder and thrust vectoring control powers at 40,000 ft and at low 
speed represent the maximum available control power. In this case, it is assumed that the vertical tail is 
not affected by the high angle of attack and that the engine is at maximum afterburner. It is concluded 
that the thrust vectoring control power is more effective, especially at the higher thrust settings, at the 
higher altitudes, and at the lower speeds where the dynamic pressure is low. 

For steady-state flight, when vectoring occurs, a decrease in forward thrust also occurs. This 
decrease is to a function of the thrust deflection angle and vectoring efficiency. The decrease is transitory 
and comparable to the drag and deceleration induced when the rudder deflects for directional control. 
As a result, for no compensation in net or gross thrust is made during the stability analysis. 

STABILITY ANALYSIS USING SMALL PERTURBATION THEORY 

Dynamic stability and control depends on the aerodynamics, mass and inertia, thrust characteristics, 
and flight conditions of the aircraft. Small linear perturbation theory was used to analyze the effects on 
the lateral-directional characteristics when thrust vectoring was incorporated with various tail sizes. 

Small perturbation theory mathematically defines the aircraft response as the change of time of the 
motion variables relative to a steady-state flight condition caused by an external or internal disturbance 
(ref. 31). The success of the small perturbation theory is caused by the fact that many aerodynamic 
forces are linear functions of disturbances, and disturbed flight of considerable violence can occur with 
small values of linear or angular disturbance velocity (ref. 32). The derivation of the linear perturbation 
equations is presented in appendix A. 
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Engine Thrust Setting Considerations 


Three control power settings were used in this study. These settings include maximum afterburner, 
military power, and steady level flight thrust. For all but the specific speeds, maximum afterburner and 
military power thrust settings would cause substantial accelerations in the forward direction and pitch 
attitude. This acceleration would imply violation of the small perturbation theory assumption of steady 
level flight. On the other hand, using the maximum afterburner control power at all flight settings would 
be representative of the maximum available control power. Small perturbation theory states that the 
total velocity is the sum of the steady state and disturbance velocities. The disturbance velocity can be 
considered the disturbance acceleration multiplied by an increment of time. If the increment of time 
becomes infinitely small, then the disturbance velocity also becomes infinitely small. Then, looking 
at the small perturbation results for maximum afterburner and military power settings at all speeds 
is representative of the maximum available control power at an instantaneous time t or an increased 
control effectiveness when compared to steady-state thrust Observing the effect of increased control 
effectiveness of postexit vanes at maximum afterburner would be indicative of the effects of an increased 
ratio of vane to rudder deflection scheduling or the increased effectiveness of an axisymmetric balanced 
beam pitch/yaw vectoring nozzle. 

The differences in thrust vectoring control power because of engine thrust setting for 40,000 ft are 
shown in figure 37. At the lower and higher Mach numbers, the control power for the steady-state flight 
engine setting is nearly that for maximum afterburner, especially at higher altitudes. The lower density at 
higher altitudes requires more thrust at slower speeds because of the increased aircraft trimmed attitude. 
For most of the frequency response plots to be presented, however, the Mach number is 0.9 at 40,000 ft 
to show the differences in magnitude and phase because of different thrust settings. 


Small Perturbation Results 


The transient response characteristics from the small perturbation analysis consist of the spiral, roll, 
and dutch roll modes and were calculated using the STAB1 program (appendix B). The denominator 
of the lateral-directional transfer function of the aircraft, the characteristic equation, provides real roots 
that are the inverse of the roll and spiral mode time constants, while the complex roots provide the 
dutch roll frequencies and damping ratios (appendix A). These open-loop modes are a function of the 
stability coefficients of the aircraft and are not dependent on control power. 


The controllability of sideslip, roll, and heading is a result of the aerodynamic or thrust vectoring 
control power and is located in the numerators in the aircraft open-loop transfer functions. The open- 
loop transfer functions are and Like the denominator of the transfer function, the 

characteristic equation, roots of the numerators provide time constants, while the complex roots provide 
the frequencies and damping ratios (refs. 31 and 32). 


When comparing the open-loop APAS results with those obtained from YF-16 aircraft flight tests, 
there was a very close correlation at Mach 0.2 for the dutch roll frequency and the roll time constant 
(refs. 26 and 27). The dutch roll damping and heading time constant did not match with those of the 
YF-16 aircraft. This difference would be expected because these numbers are usually very small and 
hard to predict. The results of the stability analysis followed the general trends when compared with 
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altitude and Mach number flight test data. As the Mach number increased, so did the differences between 
the STAB1 results and the real F-16 aircraft. Because the open-loop flying qualities are different than 
those of the real F-16 aircraft with feedback loops, this study is only a qualitative investigation of the 
trends associated with thrust vectoring and not a quantitative study. The effects on the aircraft caused by 
the additional mass and inertia, the changes in altitude and speed, the addition of thrust vectoring yaw 
power to the rudder control power, and the decrease of tail size are presented in figures 38 through 64. 

Thrust Vectoring Control System Mass and Inertia 

The addition of the inertia and mass of the thrust vectoring system for the baseline aircraft had 
an effect on the aircraft characteristics. The rolling and spiral modes were not significantly affected by 
the TVCS addition (figs. 38 and 39). The dutch roll frequency at lower speed was not affected by the 
changes in masses and inertias. As the speed increased, however, the frequency steadily decreased to 
11 percent at Mach 1.8. In addition, the dutch roll damping increased with low speed and decreased with 
higher speeds as shown in figures 40 and 41. The bank angle frequency, like the dutch roll frequency, 
was also affected by the addition of the thrust vectoring control system (fig. 42). 

Altitude 

There are large changes in the flying characteristics because of the increase in altitude for the 
characteristic equation for 100 percent tail configuration with a TVCS. The spiral mode root shows a 
very small increase in magnitude at higher Mach numbers cause by altitude that indicates a decrease in 
the time constant (fig. 43). The roll stability mode root markedly decreases with increasing altitude as 
do the dutch roll mode frequency and damping (figs. 44 to 46). The errors in the lower Mach numbers 
for increasing altitudes do not represent steady-state flight and the magnitudes of the spiral roots and 
dutch roll damping ratios are so small that they will be very inaccurate. 

The open-loop frequency response plots for the aircraft show a large decrease for the magnitude 
measured in decibels at all frequencies for the sideslip angle transfer function. For all the response 
plots, the expected shift in damped natural frequency with an increase of altitude was observed. For the 
heading, bank-angle, and open-loop transfer functions, losses in magnitude occur at frequencies above 
the damped natural frequency (figs. 47 to 52) and (ref. 33). 

Steady-Level Flight and Maximum Afterburner Thrust 

The effects of thrust vectoring at various engine settings was investigated. The larger effects of 
lateral-directional thrust vectoring are at the higher altitudes and lower speeds, where the rudder control 
power is reduced, and at higher thrust settings, where the thrust vectoring control power is largest. 
At higher altitudes, the control power level for steady level flight approaches that of the maximum 
afterburner thrust setting for high-speed flight caused by increased thrust requirements and at the very 
low speeds because of the induced drag from higher angles of attack. At the mid-subsonic Mach 
numbers, the maximum differences in vectoring power occur (fig. 37). 
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Adding the steady-state thrust vectoring control power to the rudder control power for the sideslip 
angle, open-loop transfer function has a significant magnitude increase that is comparable to the change 
of magnitude from sea level to an altitude of 40,000 ft (figs. 53 and 47). Figure 54 displays that at 
lower frequencies there was positive shift in the phase, especially for maximum afterburner and at values 
above the damped natural frequency where there was no differences in phase. 

The combination of the maximum afterburner thrust vectoring side force and yawing moment in 
addition to the rudder side force, yawing and rolling moments caused substantial decreases in the overall 
magnitude and a shift in phase for the heading angle frequency response (figs. 55 and 56). The addition 
of maximum afterburner from steady-state thrust vectoring causes the three roots of the numerator to 
go from two stable and one unstable to all three at the origin. This shift in roots to the origin is for 
maximum afterburner vectoring with rudder for configuration A and occurs throughout the entire stable 
speed regime. Because this phenomenon occurs at steady-state thrusts for high- and low-speed flight, 
careful attention must be paid to thrust vectoring for both steady-state and accelerating flight. 

The primary cause of this substantial change at lower frequencies is because of the addition of 
the thrust vectoring yaw moment. The heading angle describes the direction of the flight path and is 
affected in the positive direction by the side force and in the negative direction by the yawing moment. 
For low frequencies, yawing moment ( N ) from thrust vectoring causes large decreases in the magnitude. 
Meanwhile, the addition of thrust vectoring side force ( Y ) results in increases in the magnitude. As the 
frequencies approach the damped natural frequency the effects are opposite, the larger yawing moment 
increases magnitude, and the larger side force decreases the magnitude. The large decreases in magnitude 
caused by the relatively small addition of thrust vectoring side force and yawing moment suggests that 
careful attention should be paid to heading angle response when designing for precision tracking flight 
maneuvers at high altitudes and at high thrust situations such as close-air formation flying and in-flight 
refueling. 

The bank-angle frequency response was significantly affected with the introduction of maximum 
afterburner yaw vectoring (fig. 57). A large increase in the magnitude at lower frequencies was caused 
by both steady-state and maximum afterburner thrust vectoring. The differences between maximum 
afterburner, steady-state thrust, and rudder only controls indicate a complex relationship for the bank- 
angle response, especially at higher frequencies. The dip in the magnitude for the standard rudder and at 
steady-state thrust vectoring is caused by the open-loop transfer function acting like a aileron response 
with a complex root rather than a rudder response with no complex roots as the rudder and maximum 
afterburner thrust vectoring combination acted. This difference in the response characteristics is most 
evident in the phase response, where there is a difference of 180° at the higher frequencies (fig. 58). The 
large changes in the bank-angle response with rudder and vane deflections suggest that modifications 
are needed for the aileron/rudder interconnect to compensate for the vectoring effects. 

Reduction in Vertical Tail Size 

Aircraft designers have suggested that the installation of TVCS increases the amount of control 
power and allows a size reduction of the vertical tail. The size reduction of the vertical tail also decreases 
the lateral-directional stability of the aircraft as well as control power of the rudder thereby changing 
the dynamics of the aircraft. To estimate the effects of vertical tail size reduction with TVCS control 
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power, the dynamic response for the aircraft was studied with maximum afterburner and deflection thrust 
vectoring at 100, 75, 50, and 0 percent tail height and compared with the baseline aircraft. 

The addition of thrust vectoring dramatically increases the magnitude characteristics of the sideslip 
transfer function (figs. 59 and 53). As expected, reduction of vertical tail height shows the natural decay 
of the damped natural frequency to the point where it does not exist on the unstable configuration D 
(figs. 59 and 60). What was unexpected was the crossover between the magnitudes of configurations B 
and C, where the smaller tail height of configuration C has a larger magnitude at frequencies below the 
damped natural frequency. The results show that there are ratios of vertical tail size to thrust vectoring 
control which yield greater magnitudes at specified frequency ranges for the sideslip transfer function. 
The reduction of vertical tail size affects the lateral-directional stability and control of the aircraft, 
while the TVCS adds directional control. These differences in the vertical tail and TVCS cause the 
disproportional changes in aircraft dynamics. 

The heading angle response for the various aircraft configurations shows the expected reduction 
in the damped natural frequencies with the reduced tail size and no natural frequency for the unstable 
configuration C (figs. 61 and 62). The magnitude responses increase at lower frequencies with a re- 
duction of tail size from configuration A to configuration C which corresponds with the responses for 
the reduction of a decrease of thrust as seen in figure 55. The magnitude responses crossover above 
the damped natural frequencies. The baseline aircraft with rudder control has the largest response. The 
phase responses for the heading angle show that large changes can occur with the addition of thrust 
vectoring control power and the reduction of tail size as the root locus shift poles from real to imaginary. 
Like the sideslip responses, the heading responses suggest that there are ratios of thrust vectoring control 
power to vertical tail size that yield higher responses at specified frequencies. 

The bank-angle magnitudes and phase responses show similar characteristics to the sideslip and 
heading angle transfer functions (figs. 63 and 64). The complex roots in the denominator of the transfer 
function, which make up the damped natural frequency and damping, decrease in frequency with a 
decrease in tail height and are nonexistent in configuration D. The complex roots of the numerator that 
appear in the baseline aircraft are real roots in configuration A with thrust vectoring and are complex 
and increasing in frequency with a decrease in tail height. The response for lower frequency once 
again has crossovers and disproportional magnitudes for decreasing tail height. The large shift in phase 
response with the addition of TVCS to the baseline aircraft is caused by the change of the complex 
roots in the denominator to real roots. The complex and real roots in the numerator and denominator 
can significantly change the phase and magnitude characteristics of the bank-angle response. 

To meet the level 1 flying qualities for a highly maneuverable fighter aircraft, certain aircraft 
characteristics must be maintained. The dutch roll damping and frequency, spiral mode, and bank mode 
change with configuration changes. Such changes are especially pronounced at high Mach numbers and 
altitudes (ref. 31). The dutch roll frequency must be greater than 1 rad/sec. Because reducing the tail 
size reduces the dutch roll frequency, especially at higher altitudes and Mach numbers, sideslip angle, 
yaw, and roll rate feedback loops would need to be introduced (ref. 31). The dutch roll damping, roll, 
and spiral mode roots can also be maintained through the use of feedback loops like those currently 
employed on the F- 16 aircraft. Because this study did not employ these feedbacks, the exact magnitudes 
of the flying characteristics and the amount of tail size reduction cannot be predicted for a current flying 
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vehicle. Without feedback loops, tail height for an F-16 aircraft may be reduced to the difference 
between configurations C and D for steady-state flight and still be stable. 

SUMMARY 

Effectiveness and efficiency of the postexit vane TVCS are highly dependent on engine thrust 
setting, nozzle pressure ratio, and vane area, placement, and deflection. The large, simultaneous, two- 
vane deflections caused thrust deflections which forced the thrust vector to impinge on the inactive 
vane. In addition, the deflections caused large losses in effectiveness and efficiency. Such losses were 
especially pronounced at higher engine thrust settings and nozzle pressure ratios. In addition, mass 
and inertia increases caused by the thrust vectoring control system and counterbalance weight are much 
larger in magnitude than such increases for the removal of the vertical tail for a postexit vanes system. 

Reduction of vertical tail size decreased the skin friction drag coefficient and stability coefficients 
as predicted. On the other hand, the drag and weight increases of the postexit vanes hardware eliminated 
the potential drag and weight reductions except for the gains to be made through use of optimum pitch 
vectoring. 

Compared with postexit vanes, future thrust vectoring systems will have greatly improved efficiency, 
effectiveness, and drag and weight effects. Thrust vectoring will be a benefit of control when the 
aircraft is in the spin regime. Conversely, such vectoring will not affect the open-loop susceptibility of 
the aircraft to entering into a spin as would a reduction in the size of the vertical tail. Thrust vector 
becomes significant at higher altitudes and lower speeds, where the dynamic pressure is low. The thrust 
vectoring control power available is a function of the gross thrust and not the net thrust for all Mach 
numbers above 0.4 to 0.5 because of the effects of ram drag. 

Adding thrust vectoring to the rudder control power and changing altitude and tail size have 
significant affects on the stability and control characteristics of F-16 aircraft Meanwhile, the addition 
of the mass and inertia of the thrust vectoring control system has very small effects. The results 
presented are the open-loop transfer functions of the aircraft and do not represent the closed-loop, 
multiple feedback system currently employed by F-16 aircraft. It is not possible to accurately state the 
amount of tail size reduction. It is possible, however, to observe that the tail size of the aircraft can be 
reduced. 

The most important parameter for observing the improved yaw characteristics is the sideslip transfer 
function because this function represents the power of the directional control. When incorporating 
maximum afterburner or steady-level thrust, significant increases in magnitude can be achieved. The 
differences in the frequency responses suggest that there is an optimum combination of thrust vectoring, 
tail size, rudder control power, and vane deflection scheduling for the highest magnitude. 

The large changes in phase angles and magnitudes for heading and roll responses with relatively 
small increases in thrust vectoring control power or decreases in tail size suggest that careful atten- 
tion should be paid to these characteristics and the integration of adding on yaw vectoring to the 
rudder control. 
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CONCLUDING REMARKS 


At high altitudes, low speeds, and high engine thrust settings, thrust vectoring control can signifi- 
cantly increase aircraft maneuverability and spin control capabilities. 

The effectiveness and efficiency of postexit vane thrust vectoring control systems are highly de- 
pendent on the engine thrust setting, the nozzle pressure ratio, and the area, placement, and deflection 
of the vanes. 

The drag and weight reductions caused by the decrease in vertical tail size are not large enough 
to compensate for the increases because of the current postexit vane thrust vectoring systems. Future 
systems will have improved efficiency, effectiveness, drag, and weight effects to increase the performance 
of the aircraft 

Large disproportional changes in the magnitude and phase characteristics of the frequency response 
for the sideslip, heading, and bank angle transfer functions occur with the introduction of thrust vectoring 
and the reduction of vertical tail height. 

The vertical tail height can be reduced significantly for an aircraft by using thrust vectoring; 
however, only through careful attention to the tail size, vectoring effectiveness, as well as vane to 
rudder scheduling will the optimum control and performance be obtained. 

For future research, it is recommended that the effects of a more advanced axisymmetric pitch and 
yaw thrust vectoring nozzle design be investigated. To increase two-vane effectiveness and efficiency, 
it is recommended that all postexit vanes have coordinated symmetric and asymmetric deflections to 
eliminate jet impingement interference and to allow inactive vanes to act as aerodynamic control surfaces. 

It is necessary to investigate the effects of a feedback control system to determine the amount 
of vertical tail height reduction allowable to maintain level 1 flying qualities. Further studies on the 
effects of pitch and yaw thrust vectoring control power during nonsteady-state flight conditions and 
nonlinear regimes, such as takeoffs, landings, high angle-of-attack flight, and maneuvering flight should 
be performed. 


Dryden Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, California, August 22, 1990 
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APPENDIX A 


EQUATIONS OF MOTION 

The following sum of the forces and moments are extracted from refs. 22 and 23: 

E F x = m(U + QW - RV ) 

E F y = m(V - PW + RU) 

E F z = m(W + PV- QU ) 

EL = I XX P + QR(I ZZ - Iyy) - I yz (Q 2 - R 2 ) - I XZ (R + PQ) - I xy (Q - PR) 

EM = IyyQ + PR(I XX - I zz ) + I XZ (P 2 - R 2 ) - I yz (R - PQ) - I xy (P + QR) 

S N = I ZZ R + PQ(Iyy - Ixx) - Ixy{P 2 - Q 2 ) - Ixz{P ~ QR) ~ Iyz{Q + PR) 

Assuming that the c.g. position is on the centerline of the aircraft 

Ixy = Iyz = 0 

EL = I XX P + QR{Jzz — Iyy) — Lcz(L- T" PQ) 

EM = IyyQ + PR(I XX - I zz ) + I XZ (P 2 - R 2 ) 

EAT = I ZZ R + PQ(Iyy — I XX ) — Ixz(P — QR) 


The lateral-directional aircraft characteristic forces and moments are 

E F y = (Cy 0 p + Cy P ^j + c y r + C ys 6 ^ V s + m 9 sin <t> cos 9 

( Ph Rb \ 

+ °IP2U + Ce -2U + Cl ‘ 6 ) « S6 
EJV = (CnpP + C np ^ + C„ r ^| + C„ S S) qSb 
where 6 is the control surface deflections for ailerons, rudders, and thrust vectoring. 
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Adding the equations of aircraft motion gives 


m(V-PW + RU) = (cyeP + Cy^ + ijS C, r ^ + C ys S^qS 

+mg sin 0 cos 6 

I XX P + QR(I ti - Iyy) - I XZ (R + PQ ) = {c tf l 3 + c lp ^ + c tr i£ + C ls sj qSb 
IzzR + RQ(Iyy ~ I xx ) — Ixz(P ~ QR ) = ^ CnpP + + ^ nr 2 U^~ 


Assuming small perturbations about a reference flight condition, the velocity components can 
be written 

U = U 0 + u 

V = V 0 + v 
W = Wq + w 

P = P 0 + p 

Q = Q o + 9 

R = Rq + r 

where 

Vq = Wq = Pq = Qo = Rq = 0 

for steady-level flight. Assuming that the perturbations are small, the products would be approx- 
imately zero 

m(v + rU 0 ) = (Cypp + C Y P ^+ Cyr W^ + CysS ) S + m9 sin ^ cos 6 
IxxP ~ Ixzr = (c i0 P + C £p ^j- + C ir ^- + Q/) qSb 
Izzr - I XZ P = (Cnp/3 + C np ^ + Cn r ^ + qSb 


Assuming sin 0 « 0, dividing through by the mass for the forces and inertias for moments and 
dimensionalizing the coefficients are as follows: 

v + tUq = 50 cos 6 + Yp/3 + Y p p + Y r r + Y$5 
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where the constants 


p - Air = Lpfi + L p p + L r r + Lg8 
r - B\p = Npt 3 + N p p + N r + N$6 


A\ = 


B 1 = 


1 XZ 

Ixx 

Ixz 

Ixx 


taking the Laplace transform gives 

(sU 0 - Y p )P(s) - (sY p + g cos 0)<f>{s) + s(U 0 - Y r )*(s) = Y s 6(s) 


-LpP(s) + (s 2 - (f>(s) - ( s 2 A\ + sL r ) tf(s) = L 6 8{s ) 

- {N p + Nfr) /3(s) - (s 2 Bi + sN p ) <f>(s) + (s 2 - siV r )tf (s) = N s S(a) 
putting these equations in matrix form yields 


C sUo-Yp ) -(sYp + g cos 9) s(U 0 -Y r ) 


Pis) ) 

<5(s) 

/ 

-Lp(s) +(s 2 - L p s) ~(s 2 Ai + sLr) 

< 

3(i) 

► = < 

_-(N p + N Pt ) - (s 2 B\ + sN p ) +[s 2 -sN r ) 


J 

\ 


Solving for the transfer functions yields the following relationships: 


8{s) 


Y s -{sYp + g cos 9) +s{U Q - Y r ) 

Ls + (s 2 - LpS ) - (s 2 Ai + sLr) 

N s - ( s 2 B l + sN p ) + (s 2 - sN r ) 


sUo - Yp) -( sY p + g cos 9) s{Uq - Y r ) 

—Lp(s) +(s 2 - L p s ) - (s 2 Ai + sLrj 

- (N 0 + Npj) - (. s 2 Bi + sN p ) +(s 2 - sN r ) 


<f)(s) 


( sUo-Yp ) Y 6 +s(U 0 -Y r ) 

- Lp(s ) Lfj ~(s 2 A\ + sLr ) 

_-(Np + Np T ) N s +(s 2 — sN r ) 


8(s) 

(sUo - Yp) -( sYp + g cos 9) +s(Uq - Y r ) 

- Lp(s ) +(s 2 - Lps) -( s 2 j 4 i + sLr) 

_ - (Np + Np T ) - (s 2 Bi + sN p ) +(s 2 -sN r ) 
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(sU 0 - Yp) - ( sYp + g cos 9) Y s 
-Lp(s) +(s 2 - L p s ) Ls 

ff(g) ~ (Np + Nfo) - (s*B\ + sNpj Ns * 

6(s) ' ( sU 0 -Yp ) ~(sY p + geos 9) +s(U 0 -Y r ) 

— Lp(s ) + (g 2 — Lps) — (s 2 j4i + g£r) 

- (Np + iV^) - (s 2 #! + sATp) +(g 2 - sN r ) 

with the transient response characteristics being the denominator 

( sUo-Yp ) -(sY p + gcos9) +s(U 0 -Y r ) 

A = ~Lp{s) +(g 2 - Lps) ~(s 2 A\ + sLr) 

-(Np + Np T ) -(s 2 B x +sNp) +(s 2 — sN r ) 

Solving this matrix reveals one root indicating neutral heading stability, one first-order root known 
as the spiral mode, one highly convergent root known as the rolling mode, and one slowly damped, 
low-frequency, oscillatory, dutch roll mode. These roots are expressed as 

/3(g) _ s [ Aps 3 + Bps 2 + Cps 1 + Dp j 
J{s) ~ s[^4s 4 + £s 3 + Cs 2 + Ds 1 + E ) 

(f)(s) 5 [A<ps 2 + B^g 1 + Cj 

5(st = TiA^^Rs^C^T^^E] 

'Jr(g) g L/4<[r5^ -f- BtyS 2 + Cys 1 -|- .D\|rj 
~6{s) = ~4A? + Bs 3 + Cg 2 + D? + E] 

Solving for the roots of the equation results in the following poles and zeros: 

/3(g) kp (g + £i) (g + i 2 ) (f + h) 

6( s ) (g + Pi)(g + P 2 ) ( fi2 + 2 Cw + w 2 ) 

</>(g) _ ( g2 + 2 C^ + ^|) 

^( s ) (g + Pl)(g +P 2 ) (s 2 + + ^ 2 ) 

\P(s) (g + z\) (g 2 + 

<3( s ) (s + p\)(s + P 2 ) (g 2 + 2£uj + o> 2 ) 

Placing the results in standard aircraft format 

/3(g) _ k J_ ( + Q G§ + Q (t + 0 

* (s) (£ + 0(A + 1 )(£ + 2 + 1 ) 
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It is important to note that the standard format implies that the roll and heading responses will have 
complex roots. On the other hand, certain configurations drive all the roots of the numerators to the 
real axis, and the standard form is no longer valid. 
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APPENDIX B 

COMPUTER PROGRAMS 


Postexit Vane 

1. TVDAT-The 11 x 11 matrix data tables from the NASA Langley cold-jet tests for each vane, 
nozzle pressure ratio, and engine setting. 

2. JET-A series of six programs designed to calculate the thrust deflection, efficiency, and di- 
rection for each vane, vane deflection combination, NPR, and engine setting studied. Output is in 
EZPLOT format. 

3. MOD-A series of six programs designed to simulate the results of the JET series for deflection 
and efficiency of the one vane and simultaneous vane-deflection cases at various NPRs, engine settings, 
and vane sizes and placements. Output is in EZPLOT format. 


Aerodynamic 

1. APAS.GEOM-The file for the APAS II program where the F-16 aircraft geometry was input. 

2. F16.0UT-0utput file from the APAS run. 

3. TOSS/PILE-Two series of programs that compile aerodynamic coefficients into a series of data 
files from the APAS output file F16.0UT. Output files are in EZPLOT format. 


Engine 

TOSS5-A program that compiles the thrust and NPR results from the output of the F-100 EMD 
SPEC for both military power and maximum afterburner at various altitude and Mach numbers into a 
series of data files. Data files are in EZPLOT format. 


Stability 

1. STABl-A program that uses linear perturbation theory to calculate the flying characteristics and 
incorporates the aerodynamic coefficients from the PILE program, the engine output files from TOSS5, 
the induced thrust effects from windtunnel tests, and the MOD series of programs as subroutines. Most 
output is in EZPLOT format, while other output goes into a plotting or test format. 

2. BODE-An interactive program designed to produce the magnitude and phase characteristics 
for the frequency response. Output is formatted for use in EZPLOT. 
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FIGURES 



Figure 1. Three-view drawing of an F- 16 Fighting Falcon. 
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Figure 3. Single-vane deflections at military power and nozzle pressure ratio = 2. 
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Figure 4. Three- and four-vane thrust vectoring control system arrangements. 
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Figure 5. The F/A-18 high alpha research vehicle vane-size comparison. 
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Figure 6. Single-vane deflections at maximum afterburner and nozzle pressure ratio - 2. 
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Figure 7. Large vane deflection at military power for various nozzle pressure ratios 
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Figure 8. Simultaneous vane deflections at military power and nozzle pressure ratio = 2 
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Figure 9. Simultaneous B and C vane deflections at military power for various nozzle pressure ratios. 
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Figure 10. Simultaneous two-vane deflections at military power and nozzle pressure ratio = 6. 
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Figure 11. Simultaneous two-vane deflections at maximum afterburner and nozzle pressure ratio = 2. 
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Figure 15. Single-vane efficiencies at maximum afterburner and nozzle pressure ratio = 2. 
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Figure 16. Efficiencies for simultaneous two-vane deflections at military power and nozzle pressure 
ratio = 2. 
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Figure 17. Efficiencies for simultaneous two-vane deflections, at military power, and for various nozzle 
pressure ratios. 
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Figure 18. Aerodynamic Preliminary Analysis System II configurations of the F-16 Fighting Falcon. 
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Figure 29. Variation of side force coefficient with rudder deflection. 



Figure 30. Variation of rolling moment coefficient with rudder deflection. 
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Figure 33. Skin friction drag coefficient. 
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Figure 34. Effects of angle of attack on aircraft inertias. 






Figure 35. Dimensional yawing moments caused by the rudder and thrust vectoring control system at 
maximum afterburner and sea level. 
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Figure 36. Dimensional yawing moments caused by the rudder and thrust vectoring control system at 
maximum afterburner and 40,000 ft. 
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Figure 37. Total yawing moment control for various engine thrust settings for configuration A 
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Figure 38. Spiral mode root changes with thrust vectoring control system addition at 40,000 ft for 
configuration A. 





Figure 39. Rolling mode root changes with thrust vectoring control system addition at 40,000 ft for 
configuration A. 



Figure 40. Dutch roll frequency changes with thrust vectoring control system addition at 40,000 ft for 
configuration A. 
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Figure 41. Dutch roll damping changes with thrust vectoring control system addition at 40,000 ft for 
configuration A. 



Figure 42. Bank angle frequency changes with thrust vectoring control system addition at 40,000 ft for 
configuration A. 
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Figure 47. Altitude effects on the sideslip angle magnitude characteristics for configuration A at maxi- 
mum afterburner and M = 0.9. 
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Figure 48. Altitude effects on the sideslip angle phase characteristics for configuration A at maximum 
afterburner and M = 0.9. 
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Figure 49. Altitude effects on the heading angle magnitude characteristics for configuration A at 
maximum afterburner and M = 0.9. 
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Figure 50. Altitude effects on the heading angle phase characteristics for configuration A at maximum 
afterburner and M = 09. 
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Figure 51. Altitude effects on the roll angle magnitude characteristics for configuration A at maximum 
afterburner and M = 0.9. 
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Figure 52. Altitude effects on the roll angle phase characteristics for configuration A at maximum 
afterburner and M = 0.9. 
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Figure 53. Effects of thrust vectoring control power on sideslip angle magnitude characteristics for 
various engine settings for configuration A at 40,000 ft and M = 0.9. 
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Figure 54. Effects of thrust vectoring control power on sideslip angle phase characteristics for various 
engine settings for configuration A at 40,000 ft and M = 0.9. 
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Figure 55. Effects of thrust vectoring control power on heading angle magnitude characteristics for 
various engine settings for configuration A at 40,000 ft and M = 0.9. 



Figure 56. Effects of thrust vectoring control power on heading angle phase characteristics for various 
engine settings for configuration A at 40,000 ft and M = 0.9. 
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Figure 57. Effects of thrust vectoring control power on bank angle magnitude characteristics for various 
engine settings for configuration A at 40,000 ft and M = 0.9. 
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Figure 58. Effects of thrust vectoring control power on bank angle phase characteristics for various 
engine settings for configuration A at 40,000 ft and M = 0.9. 
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Figure 59. Effects of tail size on sideslip angle magnitude characteristics at maximum afterburner, 
40,000 ft, and M = 0.9 with thrust vectoring. 
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Figure 60. Effects of tail size on sideslip angle phase characteristics at maximum afterburner, 40,000 
ft, and M = 0.9 with thrust vectoring. 
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Figure 61. Effects of tail size on heading angle magnitude characteristics at maximum afterburner, 
40,000 ft, and M = 0.9 with thrust vectoring. 
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Figure 62. Effects of tail size on heading angle phase characteristics at maximum afterburner, 40,000 
ft, and M = 0.9 with thrust vectoring. 
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Figure 63. Effects of tail size on bank angle magnitude characteristics at maximum afterburner, 40,000 
ft, and M = 0.9 with thrust vectoring. 
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Figure 64. Effects of tail size on bank angle phase characteristics at maximum afterburner, 40,000 ft, 
and M = 0.9 with thrust vectoring. 
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